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Popular summary 

Satellite remote sensing is such a complex task that could be done only by assuming that 
there are no radiative interactions between areas that have different cloud properties. It is 
true that some recently proposed methods consider such interactions, but these novel 
methods are not yet ready for operational use. As a result, researchers must continue 
relying on one-dimensional (ID) radiative theory, which does not include the 
interactions. The hope is that the neglected three-dimensional (3D) interactions do not 
matter much in remote sensing applications, and so current methods give accurate results. 

This study addresses the question whether ID radiative transfer theory describes well 
how the clouds’ solar reflection depends on viewing angle. The statistical analysis of a 
large set of MODIS observations indicates that in oblique backward scattering directions, 
cloud reflection is stronger than ID theory would predict. After considering a variety of 
possible causes, the paper concludes that the most likely reason for the increase lies in 3D 
radiative interactions. The results’ main implication is that cloud optical depths retrieved 
at back scattering view angles larger than about 50° tend to be overestimated and should 


be used only with great caution. 
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Abstract 


This study addresses the question whether ID radiative transfer theory describes 
well the angular distribution of shortwave cloud reflection. The statistical analysis of a 
large set of MODIS observations indicates that in oblique backward scattering directions, 
cloud reflection is stronger than ID theory would predict. After considering a variety of 
possible causes, the paper concludes that the most likely reason for the increase lies in 3D 
radiative interactions. The results’ main implication is that cloud optical depths retrieved 
at back scattering view angles larger than about 50° tend to be overestimated and should 
be used only with great caution. 
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1. Introduction 


Satellite remote sensing is such a complex task that could be done only by 
assuming that there are no radiative interactions between areas that have different cloud 
properties. It is true that some recently proposed methods (e.g., Marshak et al. 1998, 
Oreopoulos et al. 2000, Faure et al. 2001, Vamai and Marshak 2002) consider such 
interactions, but these novel methods are not yet ready for operational use. As a result, 
researchers must continue relying on one-dimensional (ID) radiative theory, which does 
not include the interactions. The hope is that the neglected three-dimensional (3D) 
interactions do not matter much in remote sensing applications, and so current methods 
give accurate results. 

In recent years, several observational studies examined the validity of this hope. 
These studies found that, under certain conditions, 3D effects cause significant problems. 
Specifically, they revealed that 3D effects can make clouds appear too smooth (Marshak 
et al. 1995, Davis et al. 1997), too bright and thick (Loeb and Davies 1996, Loeb and 
Coakley 1998), and artificially asymmetric (Vamai and Marshak 2002). 

While the papers above focused mainly on overhead satellite views, some studies 
examined 3D effects for oblique views. A comparison of GOES and Meteosat radiances 
for scenes that were viewed from different directions by the two satellites did not reveal 
any influence of 3D effects (Rossow 1989). Examining ERBE (Earth Radiation Budget 
Experiment), AVHRR (Advanced Very High-Resolution Radiometer), and POLDER 
(Polarization and Directionality of the Earth’s Reflectances) data, some other studies 
found that for low sun, 3D interactions such as shadowing make clouds appear too dark 
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from oblique views facing the sun, and that this makes ID retrievals underestimate cloud 
optical thickness (Loeb and Davies 1997; Loeb and Coakley 1998; Buriez et al. 2001). 
Theoretical studies (e.g., Davies 1984; Kobayashi 1993; Vdmai 2000; Iwabuchi and 
Hayasaka 2002) have long suggested that 3D effects have an opposite influence for 
oblique views facing away from the sun— but the observations cited above have not 
confirmed unambiguously the existence of this enhanced backscatter from sunlit slopes. 
Most recently, Zuidema et al. (2003) found that in highly heterogeneous cumulus 
congestus clouds, oblique backscatter reflectances observed by MISR (Multi-angle 
Imaging SpectroRadiometer) exceeded 3D radiative transfer calculations based on cloud 
structure retrieved from the MISR nadir camera using ID theory. 

The goal of this paper is to analyze the view angle dependence of cloud reflection 
and its effect on cloud optical depth retrieval. Special focus is put on the enhancement in 
backward scattering directions. For this, the study uses high resolution MODerate- 
resolution Imaging Spectroradiometer (MODIS) observations of a wide variety of cloud 
types. First, Sections 2 and 3 describe the data we analyzed and the methodology we 
used, then Section 4 discusses the basic results. Next, Section 5 examines whether the 
observed features are really caused by 3D effects or perhaps by some other factors. 
Finally, Section 6 offers a brief summary and discusses the results’ main implications. 

2. Observations 

This study took advantage of the unprecedented abundance of high-quality, easy- 
to-use, and freely available cloud products from new Earth Observing System (EOS) 
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satellites. In particular, it used observations by MODIS on board the Terra and Aqua 
satellites. The analysis focused on the operational 1 km-resolution cloud optical thickness 
product (Platnick et al. 2003), but it also considered the calibrated reflectances at 0.86 }i m 
wavelength, the cloud phase product, and some sun-view and geolocation parameters. All 
products were generated by the Collection 4 version of the operational MODIS data 
processing software. The analysis used data from each Friday in November 2002: 
November 1, 8, 15, 22, and 29. The one- week separation between subsequent days 
reduced potential sampling biases, because it is unlikely that the same cloud systems 
could persist through several Fridays. 

For each Friday we took 10 MODIS granules from the Terra and Aqua satellites, 
which yielded a total of 5 x (10 + 10) = 100 granules. (Each granule is an image covering 
an approximately 2000 km by 2000 km area.) The 10 granules were chosen so that they 
form a ring around the Earth at roughly 40° North latitude. Terra observed this ring in the 
morning, and Aqua, in the afternoon. The histogram of solar zenith angle is shown in 
Figure 1. 

3. Methodology 

Simple schematic arguments (Fig. 2) illustrate that 3D effects are expected to 
enhance cloud reflection into backscatter directions. Since ID data processing algorithms 
do not account for this enhancement, they assume that clouds viewed from backscatter 
directions appear so bright because they are very thick. Therefore our basic approach to 
detecting 3D effects was to examine the question: Do ID retrievals yield systematically 
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larger optical thicknesses when clouds are viewed from backscatter directions? 
Obviously, the clouds’ true optical thickness (x) does not depend on the direction form 
which a satellite views them; therefore, a statistically significant increase in x can indicate 
the presence of 3D effects. (Conceivably, factors other than 3D effects could also cause 
such increases in the observations, and so Section 5 will examine whether these other 
factors also play an important role.) 

The geometry of MODIS observations implies that clouds are viewed from 
backscatter directions at one edge of MODIS images, and from forward scattering 
directions at the other edge (Figure 3). As a result, examining the view angle-dependence 
of retrieved x values is equivalent to examining how x varies across the satellite track. Let 
us note, however, that the Sun-synchronous orbits of the Terra and Aqua satellites imply 
that at 40° latitude, clouds are not viewed from the exact forward and backward 
directions, but approximately 50° off the plane of solar azimuth. 

4. Results 

Before examining the x-retrieval results, let us briefly pause at the observed 
reflectances. Figure 4 reveals that their basic behavior is qualitatively consistent with ID 
radiative transfer theory: Reflectances are lowest at nadir and increase in oblique 
directions, especially for forward scattering. However, there are also some qualitative 
differences that are consistent with 3D effects, which increase backscatter and reduce 
forward reflection. Most notably, the observed reflectances do not level off at very 
oblique view angles as ID theory predicts (Fig. 4). Unfortunately, a quantitative 
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comparison of observed and simulated reflectances would be a complicated task, because 
observational, environmental, and cloud parameters vary substantially across the 
observed scenes. Therefore we take advantage of the MODIS operational retrievals that 
already considered all these variations, and we turn our attention to the retrieved x-values. 

Figure 5 displays the view angle-dependence (i.e., cross-track variability) of the 
retrieved optical thicknesses’ mean and standard deviation. The figure was obtained by 
combining all those pixels from the 100 MODIS granules that contained liquid phase 
clouds with x > 2. Pixels with x < 2 were excluded because the number of such pixels 
varies significantly with view direction due to a simple reason that has nothing to do with 
3D radiative effects: The cloud masking algorithm (Ackerman et al. 1998; Platnick et al. 
2003) is more effective in detecting thin clouds at oblique views, because the viewing 
path through the clouds is longer (than at overhead views). 

The figure reveals significant increases at oblique backscatter directions, which is 
fully consistent with the expected influence of 3D radiative effects. We note that a similar 
figure for ice clouds (not shown) displays this trend even stronger— but because it would 
be difficult to separate the influence of 3D effects from the influence of uncertainties in 
ice crystal phase functions (Yang and Liou 1996), we limit the current study to liquid 
water clouds. 

Figure 6 indicates that a major contribution to the increase in backscatter 
directions comes from an increase in the ratio of saturated pixels. Saturation occurs when 
the observed reflectances exceed the range in which ID reflectances are sensitive enough 
to optical thickness for allowing meaningful retrievals. In such cases the retrievals do not 
return unrealistically high x values, but instead report the maximum allowed x value in 
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the range of 98-100. The increase in the number of saturated pixels is easy to understand 
in terms of 3D effects: as 3D effects make clouds brighter in backscatter views, the 
saturation threshold will be exceeded more frequently. 3D effects increasing the rate of 
saturation not only increase the mean optical thickness but also alter the structure of 
cloud fields. This implies that one needs extra caution if oblique observations are 
included into studies of cloud structure using 2-point statistics such as structure functions, 
autocorrelation functions, power spectra, fractal dimensions or intermittency (e.g., 
Oreopoulos et al., 2000). Still, Figure 6b indicates that the increase in backscatter views 
is significant (up to 30% for © > 60°) even if only non-saturated pixels (with x < 98) are 
considered. 

While the increases in backscatter directions (especially for 0 > 50°) can certainly 
come from 3D radiative effects, it is unclear whether the variations at less oblique views 
and in forward scattering views are related to 3D effects. As mentioned in the 
introduction, theoretical simulations and earlier observations showed that 3D effects tend 
to reduce reflection in forward directions, whereas Figures 5 and 6 show some increases 
with 0. We are not yet ready to correctly interpret these increases; thus the present paper 
focuses only on interpreting the sharp increases observed in backscatter directions. 

5. Influence of factors other than 3D radiative effects 

In order to make sure that the x-increases in backscatter directions are indeed 
caused by 3D radiative effects, this section examines a variety of factors that could 
conceivably also cause similar increases. Because the increase appears most pronounced 
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in Figure 6a, this section uses the ratio of saturated pixels to explore possible alternative 
explanations for the x-increase in backscatter directions. 

Surface reflection: Figure 7 indicates that the x-increase is present over both land 
and ocean surfaces. Since it is unlikely that the cloud retrievals would have such similar 
errors in the reflection characteristics of these very different surface types, one can 
exclude uncertainties in surface properties from the causes of x-increase. (We note that 
the observations do not include either sun glints or direct back-reflection hotspots.) 

Atmospheric correction: Errors in the retrieval algorithm’s handling of gaseous 
absorption and Rayleigh scattering are not likely to contribute much to the x-increase, 
since such errors would have similar effects for forward and back oblique views. 
Uncertainties in aerosol properties are also cannot be the main factor, because Figure 7 
revealed that the x-increase is strong over both land and ocean despite the typically very 
different aerosol properties over these surfaces. 

Droplet scattering phase function: Uncertainties in the phase function are unlikely 
to cause the observed x-increases, for three reasons. First, the ratio of saturated pixels (or 
the mean x value) plotted as a function of the scattering angle (not shown) does not 
indicate any increase at scattering angles around 125°, which dominate at the regions of 
the x-increase. Second, even though the Earth’s curvature and orbital parameters result in 
similar scattering angles around 0 = 30° than around 0 = 60°, the observations do not 
show any x-increases around 0 = 30°. Third, a comparison of the angular distribution of 
observed liquid and ice cloud reflectances revealed that if some ice clouds were 
mistakenly identified as liquid clouds, the mistake would act to reduce (and not enhance) 
the x-increase. 
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Local time: Since the x-increases are present in both the morning observations of 
Terra and in the afternoon observations of Aqua (Figure 8), it is unlikely that systematic 
cross-track changes in local time would play a major role in the observed x-increases. 

Latitude: Because of the orbital inclination of the Terra and Aqua satellites, 
MODIS observations have a systematic cross-track dependence: The areas observed at 
backscatter directions lie at slightly higher latitudes than the areas observed from 
overhead or forward scattering directions. To check whether this may cause the observed 
x-increases, the importance of various latitudes (A) across the MODIS track was 
equalized using the equation 



0 ) 


where N x is the number of l°-wide latitude bands in the 30°-45° range. (In this range, 
there are plenty observations at all view angles.) Figure 9 shows the x-increase in 
backscatter directions even after equalizing the importance of various latitudes across the 
track, which means that the x-increase does not arise because of a systematic latitude- 
dependence in cloud properties. 

Solar zenith angle: Similarly to the latitude, the solar zenith angle «9 0 ) also varies 
systematically across MODIS images. However, the results (not shown) indicate that the 
x-increase in backscatter directions remained strong even after equalizing the importance 
of various solar zenith angles similarly to Eq. (1). This indicates that the x-increase is a 
manifestation of 3D radiative effects that depend on the viewing and not the solar zenith 
angle. This is different from the 3D effects described in some earlier studies (Loeb and 
Davies 1996; Loeb et al. 1997; Loeb and Coakley 1998). 
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Sampling noise: Some of the variations in Figures 5 and 6 are certainly caused by 
sampling noise— that is, by the random influence of individual cloud systems. However, 
it is very unlikely that the x-increase in backscatter direction would be caused by 
sampling effects First of all, its magnitude is much larger than that of random variations 
at other view angles. Also , Figure 8 revealed that the x-increase is present in both Terra 
and Aqua observations. This is important, because Terra and Aqua observed the same 
weather systems (with only a few hours difference), but the Terra and Aqua granules 
were not aligned. As a result, a weather system observed by Terra from oblique 
backscatter directions was usually observed by Aqua from different directions. This 
implies that if an individual cloud system caused the backscatter x-increase in Terra 
observations, it would have caused a similar increase at different view angles in the Aqua 
observations. 

6. Summary and conclusions 

This study examined whether one can observe in real clouds a tendency that3D 
radiative interactions enhance cloud reflection into backscatter directions. In particular, it 
analyzed the view angle-dependence of cloud optical thickness values that were retrieved 
from MODIS observations. The idea was that if 3D effects do enhance back reflection, 
the ID retrievals (not accounting for this effect) would yield larger optical thickness 
values when clouds are observed from back scatter, rather than overhead or forward 
scatter view directions. 
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The paper used the 1 km-resolution MODIS cloud optical thickness product. The 
statistical analysis of 100 approximately (2000 km) 2 mid-latitude images from November 
2002 revealed that indeed, ID retrievals yield systematically higher optical depths when 
clouds are viewed from back scattering directions. Both the mean and the standard 
deviation of cloud optical depth distributions increased by up to 30%, and the number of 
pixels too bright for accurate retrievals jumped by up to a factor of 5 for oblique 
backscatter views. 

The results have some important implications for satellite remote sensing: 

• When the sun is less than 30°-40° above the horizon, one should exercise great 
caution in using ID cloud property retrievals based on oblique backscatter 
observations. In some studies it may be prudent to limit the use of ID retrievals to 
viewing zenith angles less than 50°. This implication can be especially important 
for geostationary satellites, which always observe some areas at such oblique 
angles. For backward viewing angles larger than 50°, one can expect a substantial 
overestimate in ID cloud optical depth retrievals. The overestimate grows 
sharply with viewing zenith angle. 

• Because 3D radiative effects greatly modify the structure of the retrieved cloud 
optical depth fields (e.g., they increase the number of very thick pixels), one 
should be very careful about including such observations into studies of 
horizontal variability, such as intermittency or fractal properties. 
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Figure captions 


Figure 1. Histogram of the used observations’ solar zenith angle. 

Figure 2. 3D effects that enhance cloud reflection into backscatter directions. The solid 
arrow indicates the paths radiation follows in a 3D scene, while dashed lines indicate the 
path radiation would follow in a ID world, (a) Trapping of radiation in 3D clouds. This 
mechanism redirects photons that— because of the droplet phase function’s forward 
scattering peak— would create intense forward reflection in ID clouds, (b) Escape of 
backscatter radiation. This mechanism enhances back reflection from sunlit slopes. 

Figure 3. A schematic view of MODIS observation geometry. 

Figure 4. A qualitative comparison of the view angle-dependence the observed 0.86 jim 
reflectances to ID simulation results. In order to minimize the effects of underlying 
surface, only the brighter 50% of observed cloudy pixels are included at each view angle. 
The ID results come from DISORT (Stamnes et al., 1988) calculations for conditions 
similar to the observations. The simulation curve is plotted above the observational one 
only for schematic comparisons. 
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Figure 5. View angle-dependence of the mean and standard deviation of x- values in the 
operational MODIS cloud product. These statistics were calculated using all liquid clouds 
with x > 2 in 100 MODIS granules. Over 60,000 cloudy pixels were used in each 
(approximately 0.1°-wide) view angle bin corresponding to a single column in MODIS 
images. 


N (0) 

Figure 6. (a) Ratio of pixels with x > 98. The ratio (R) is calculated as R(0) = — — r , 

W I>2 (©) 

where N is the number of pixels and 0 is the viewing zenith angle, (b) Mean optical 
thickness of pixels with 2 < x < 98. 

Figure 7. Ratio of saturated pixels (x > 98) over land and ocean. The ratio is calculated 
the same way as in Figure 4a. 

Figure 8. Ratio of saturated pixels (r > 98) in MODIS observations from the Terra and 
Aqua satellites. The ratio is calculated the same way as in Figure 4a. 

Figure 9. Ratio of saturated pixels (r> 98) after the influence of various latitudes across 
the MODIS track is equalized using Eq. (1). 
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